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The use of Heavy Quark Symmetry to study bottom and charmed baryons leads to important 
simplifications of the non-relativistic three body problem, which turns out to be easily solved by a 
simple variational ansatz. Our simple scheme reproduces previous results (baryon masses, charge and 
mass radii, ■ ■ ■) obtained by solving the Faddeev equations with simple non-relativistic quark-quark 
potentials, adjusted to the light and heavy-light meson spectra. Wave functions, parameterized in 
a simple manner, are also given and thus they can be easily used to compute further observables. 
Our method has been also used to find the predictions for strangeness-less baryons of the SU(2) 
chirally inspired quark-quark interactions. We find that the one pion exchange term of the chirally 
inspired interactions leads to relative changes of the At, and A c binding energies as large as 90%. 
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I. INTRODUCTION 



Prior, and at around, the time the charmed baryons existence was confirmed [Da great theoretical activitydeveloped 
in order to understand and predict their masses 0-0-tJ- magnetic moments^] ; or decay properties [(J III HI- This 
activity continued in the following decades. The discovery of the Af, baryon at CERN [9(, the discovery of most of 
the charmed baryons of the SU(3) multiplet on the second level of the SU(4) lowest 20-plet [lfjj and the claims of 
indirect evidences for the semileptonic decays of Ab and E& renewed the interest in the spectroscopy and weak 
decays of heavy baryons. There are eight lowest-lying baryons containing one heavy and two light quarks (up, down, 
or strange). The quantum numbers of the charmed and bottom baryons are listed in Table [I] 

Heavy Quark Symmetry (HQS) has proved to be a useful tool to understand the bottom and charm physics [l3l 
ITiL Il5j|. and it has been extensively used to describe the dynamics of systems containing a heavy quark c or b. For 
instance, all lattice QCD simulations rely on HQS to describe bottom systems |l6|. HQS is an approximate SU (Nf) 
symmetry of QCD , being Np the number of heavy flavors: c, b, • • •, which appear in systems containing heavy 
quarks with masses that are much larger than the typical quantities (q — Kqcd, fn u ,nid,rn s • • • ) which set up the 
energy scale of the dynamics of the remaining degrees of freedom. HQS has some resemblances, in atomic physics, to 
the approximate independence of the electron properties from the nuclear spin and mass, for a fixed nuclear charge. 
Up to corrections of the order 1 0(q/m,Q), HQS guaranties that the heavy baryon light degrees of freedom quantum 
numbers, compiled in Table HJ are always well defined. The symmetry also predicts that the pair of baryons £ and 
S* or the pair H' and S* or the pair f2 and Q* become degenerated for an infinitely massive heavy quark, since both 
baryons have the same cloud of light degrees of freedom. 

However, HQS has not been systematically employed in the context of non-relativistic quark models. Non-relativistic 
quark models, based upon simple quark-quark potentials, partially inspired by Quantum Chromodynamics (QCD), 
lead to reasonably good descriptions of hadrons as bound states of constituent quarks and also of the general features 
of the baryon-baryon interaction p"?|-|l9j|. Most of the quark-quark interactions include a term with a shape and 
a color structure determined from the One Gluon Exchange (OGE) contribution |20j and a confinement potential. 
The force which confines the quarks is still not well understood, although it is assumed to come from long-range 
nonperturbative features of QCD [2l|. Those terms do not incorporate another important feature of QCD: Chiral 
Symmetry (CS) and its pattern of spontaneous breaking. Quark-quark interaction terms derived from Spontaneous 
Chiral Symmetry Breaking (SCSB) have been taken into account for the description of the nucleon-nucleon system 
and/or the light baryon spectra pi I2I I2I I2ll2r| 

In this work, we develop a variational approach for the solution of the non-relativistic three-body problem in 
baryons with a heavy quark. Thanks to HQS, the method we propose turns out to be quite simple, and leads 
to simple and manageable wave functions 2 . We consider several simple phenomenological quark-quark interactions 
which free parameters have been adjusted in the meson sector and are then free of three-body ambiguities. 



to simp, 
H3,|2g 



1 The quantities q and mq are a typical energy scale relevant for the light degrees of freedom and the mass of the heavy quark, respectively. 

2 We use a family of wave functions for which the light degrees of freedom quantum numbers are well denned. 
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TABLE I: Summary of the quantum numbers, experimental and lattice QCD masses of the baryons containing a single heavy 
quark. 7, and SJ ght are the isospin, and the spin parity of the light degrees of freedom and S, J p are stra ngen ess and the spin 
parity of the baryon (/ denotes a light quark of flavor u or d). Experimental masses are taken from Ref. [lpj and the isospin 
averaged masses are quoted, with errors counting for the mass differences between the members of the same isomultiplet. Errors 
on the lattice QCD masses have been obtained by adding in quadratures the statistical and systematic errors given in Ref. |T^|. 



With those potentials we compute the ground states of some heavy flavor hyperons and show how our simple approach 
reproduces the results obtained in Ref. [2^| by solving involved Faddeev type equations. Besides, we also study the 
£*, S', S* and Q* baryon ground states, which were not considered in the work of Ref. [28j . After that, and within 
our variational framework, we work out the predictions for the strangeness-less charmed and bottom baryons of the 
SU(2) chirally inspired quark-quark interaction of Ref. [23, applied with great success to the meson sector in Ref. [29|. 
and study the effects in these systems of including a pattern of SCSB. Some preliminary results have been presented 
in Ref. Hfl. 

The paper is organized as follows. After this introduction, we study the hamiltonian of the system, different inter- 
quark interactions (Sect.^) and the general structure of the variational baryon wave-function fSect. lllT)) . Some static 
properties (mass and charge densities, • • •) are revisited in the Sect. HVl and finally in the Sects. IVl and I VII we present 
our results and compile the main conclusions of this work, respectively. 



II. THREE BODY PROBLEM 



A. Intrinsic Hamiltonian 




In the Laboratory (LAB) frame (see Fig. QJ, the Hamiltonian (H) of the three quark (q,q',Q, with q,q' = I or s 
and Q — c or b) system reads: 

/ T72 \ 

(1) 

where m qi m q i and vtlq are the quark masses, and the quark-quark interaction terms, Vy, depend on the quark spin- 
flavor quantum numbers and the quark coordinates (x\, and Xh for the q, q' and Q quarks respectively). The nabla 
operators in the kinetic energy stand for derivatives with respect to the spatial variables xi,X2 and Xh- To separate 
the Center of Mass (CM) free motion, we go to the heavy quark frame (R, fi, r 2 ), 

ft m q xi + m q >x 2 + mQXh 

H = 

m q + m q i + niQ 

n = xi-$ h 

?2 = x 2 -x h (2) 



3 



where R and fi (fa) are the CM position in the LAB frame and the relative position of the quark q (q') with respect 
to the heavy quark Q. The hamiltonian now reads 

H = -^S + Hint (3) 



2M 

V? _ V| _ Vi ■ V ; 
2/xi 2/^2 mQ 



+VQ g (fi, spin) + VQq>(f2, spin) + £ m, (4) 

where M = (m q + m q i +mq), /x 1-2 = (l/m qq i +l/mg) 1 and Via — d/df ly ? 2 . The intrinsic hamiltonian iJ lnt 
describes the dynamics of the baryon and we will use a variational approach to solve it. It can be rewritten as the sum 
of two single particle hamiltonians {h s i p ) 1 which describe the dynamics of the light quarks in the mean field created 
by the heavy quark, plus the light-light interaction term, which includes the Hughes-Eckart term (Vi • V 2 )- 

H int = £ hr + V q An-f2,spin)-^^+ m * ( 5 ) 

i=q,q' Q i=q,q',Q 

v 2 

K = -57- + V Qi {fi,spin), i = q,q (6) 



B. Quark— Quark Interactions 

• Phenomenological Interactions: We compile here some phenomenological quark-antiquark potentials fitted to a 
large sample of meson states in every flavor sector. The general structure is as follows ( i, j = I, s,c,b): 



2 iJ2 




K (l — e -r/rA [ p-r/rg cy p -r 2 /x 

V${r) = f 1 + Xr p - A + < ao ^— e -^- + -£-«' (l - e~^A ^ 

with a the spin Pauli matrices, rrii the quark masses and 



OjOj ( 7 ) 



. . . / 2m i m i \ B . . 

\rrii + m j / 

We have examined five different interactions, one suggested by Bhaduri and collaborators [2]} (BD) and four 
suggested by B. Silvestre-Brac and C. Semay OSS (AL1, AL2, API y AP2). The parameters of the different 




FIG. 1: Definition of different coordinates used through this work. 
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TABLE II: Parameters, from Refs. |27ll28|. of different phenomenological potentials (Eqs. J7J and 10). 



potentials are compiled in Table ITTI The potentials considered differ in the form factors used for the hyperfine 
terms, the power of the confining term (p = 1, as suggested by lattice QCD calculations 32], or p — 2/3 which 
for mesons gives the correct asymptotic Regge trajectories 33]), or the use of a form factor in the OGE Coulomb 
potential. 

The usual V qq — V? 9 /2 prescription, coming from a XiXj color dependence (A are the Gell-Mann matrices) of 
the whole potential, has been assumed to obtain the quark-quark interactions from those given in Eq. ■ 

All those interactions were used in Ref (2^] to obtain, in a Faddeev calculation, the spectrum and static properties 
of heavy baryons. Our purpose is to show that our simpler variational method gives equally good results for all 
observables analyzed in [23 and, besides, it provides us with manageable wave functions that can be used for 
the evaluation of further observables. 



• Chiral Quark Cluster Interactions: Working in the context of the SU(2) linear sigma model, Fernandez and 
collaborators have derived in Ref. [23] a quark-quark interaction that, apart from the effective OGE and confining 
terms, contains a pseudoscalar (V PS ) and a scalar (V s ) potentials provided by the exchange of Goldstone bosons. 
This model was used in Ref. [2£| to analyze the strangeness-less meson spectrum with good results. We analyze 
this potential in order to check the predictions of the SU(2) chirally inspired quark-quark interactions in the 
strangeness-less heavy baryon sector. The study of baryons with strangeness would require to extend the model 
of Ref. [2^] to three flavors. A linear realization of SCSB will involve, not only the exchange of new Goldstone 
bosons, as kaon or eta mesons, but also more than one scalar meson. We are not aware of any published potential 
with such characteristics and that performs well in the meson sector. Therefore, we prefer to work in this latter 
case with baryons without strange quark content. 

In the quark-quark sector for u and d flavors, the potential reads: 
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(9) 
(10) 

(11) 
(12) 

(13) 



with t the isospin Pauli matrices, Y(x) = e x /x and ^A^Aj^ takes the value —8/3 for qq pairs in a baryon. 

The parameters given in Refs. [H Hj are: r = 0.145 fm, a s = 0.7, A CSB = 3.15 fm-\ a c = 140 MeV/fm, 
m u = md — 313 MeV, a c h — 0.027569, m^ = 138 MeV and m a = 675 MeV. To reproduce the experimental 
p — 7T mass splitting 771 — 138 = 633 MeV), the parameter r needs to be slightly reduced to 0.1419 fm. 



5 



Besides, the overall origin of energies, a&, takes the value of 122.1 MeV. Note that yXiXjj takes the value —16/3 

for qq pairs in a meson and that the isospin structure {nfj) should be multiplied by a factor —1 in the qq case. 

To use the above interaction to predict strangeness-less heavy baryon masses, such potential needs to be sup- 
plemented by interactions between heavy quarks (c and b) and light quarks (u and d). If we look at the 
phenomenological potentials discussed in Eq. Q and Tabled we see the AL1 interaction uses a light quark 
mass very close to that used in Refs. [23, (315 versus 313 MeV) and furthermore in both cases a linear 
confining potential is used. Besides, the AL1 potential provides masses for the D,D*,B,B* quite close to the 
experimental values, and it and p masses (138 and 770 MeV respectively) in good agreement with the values 
obtained from the interaction of Eq. J^J- Thus, we will compute the masses of the A c f,,E Cj 6 and S* b baryons 
with a new potential, which we call ALlx, defined as follows: 

— The light-light quark (u and d) interactions will be described by the potential of Eq. (JSJ , which incorporates 
a pattern of SCSB, and the parameters given below Eq. I|13|l. but with m u = rrid = 315 MeV, ro = 0.1414 
fm and a& = 122.4 MeV. These two latter parameters have been very slightly modified to obtain 138 and 
770 MeV (values provided by the AL1 potential) for the ir and p masses respectively, when a u, d quark 
mass of 315 MeV is used. The new set of parameters provides states with masses of 138, 1236 and 1918 
(770, 1558 and 2173) MeV with n (p) quantum numbers, in perfect agreement with the results of Ref. [29| . 

— To describe the interaction between the heavy quarks and the u and d quarks we will adopt the AL1 model. 

In this way, and by comparing results obtained from both the AL1 and AL1% interactions, we will be able to 
test the effects of the inclusion of a SCSB pattern to describe the light-light quark interaction in the u, d sector. 
Such a study, though it has already been performed in other systems, as light meson or baryon spectrum, etc, 
has not been ever carried out for heavy baryons. It is of great interest since the quantum numbers of the light 
degrees of freedom in a A— or E— type heavy baryon have a clear correspondence to those of the ir or p meson. 

III. VARIATIONAL WAVE FUNCTIONS 

In a baryon, the singlet color wave function is completely anti-symmetric under the exchange of any of the three 
quarks. Within the SU(3) quark model, we assume a complete symmetry of the wave function under the exchange 
of the two light quarks (u, d, s) flavor, spin and space degrees of freedom. On the other hand, for the interactions 
described in the previous subsection, we have that both the total spin of the baryon, Sb = (<J q + o q < + <Jq) /2, and 
the orbital angular momentum of the light quarks with respect to Q, L, defined as 

L = lx + h, with l k = —i f k x V fe , k = 1,2 (14) 

commute with the intrinsic hamiltonian. We will assume that the ground states of the baryons are in s-wave, L = 0. 
This implies that the spatial wave function can only depend on the relative distances r\, r-i and r\i = \f\ — ri\. 
Note that when the heavy quark mass is infinity (rriQ — > oo), the total spin of the light degrees of freedom, S^ght = 
(<7 q + ffq>) /2, commutes with the intrinsic hamiltonian, since the gq -o q . q i /(mQm q ^ q i) terms vanish in this limit. With 
all these ingredients, and taking into account the quantum numbers of the light degrees of freedom for each baryon, 
compiled in Tableland that in general are always well defined in the static limit mentioned above, we have used the 
following wave functions in our variational approach 3 

• A— type baryons: 1 = 0, Slight = 

\A Q ;J=±,Mj) = {|00) / ®|00) Sllght }< Q (r 1 ,r 2 ,r 12 )® \Q;Mj) (15) 

where ^fi Q (t"i, ra, n.2) = ^fi Q (?"2, fi, ^12) to guarantee a complete symmetry of the wave function under the 
exchange of the two light quarks (u, d) flavor, spin and space degrees of freedom, and finally Mj is the baryon 



3 An obvious notation has been used for the isospin-flavor (\I,Mj)j, \ls) or \sl)) and spin (\S, A^s)si ig ht ) wave functions of the light 
degrees of freedom. 
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total angular momentum third component. Note, that SU(3) flavor symmetry (SU(2), in the case of the A<; 
baryon) would also allow for a component in the wave function of the type 



E (^\M Q MsMj){\00)j® \lMs)s Usht }eu Q (ri,r 2 ,n 2 )® \Q;M Q ) (16) 

M s Mq 



with Q U Q (n, r 2 , riz) — — Qu Q (r 2 , fi, i"i 2 ) (f° r instance terms of the type r\ — r 2 ) and, the real numbers 
(jijijlmiTnam) = {jimij 2 m 2 \jm) are Clebsh-Gordan coefficients. This component is forbidden by HQS in 
the limit tuq — > 00, where Slight turns out to be well defined and set to zero for Ag— type baryons. The most 
general SU(2) Aq wave function will involve a linear combination of the two components, given in Eqs. I|15|l and 
Neglecting 0(q/iriQ), HQS imposes an additional constraint, which justifies the use of a wave function of 
the type of that given in Eq. (|15|) with the obvious simplification of the three body problem. 

£ and Y<* — type baryons: 1 = 1, Slight = 1- 

|S Q ; J = i, Mj; M T ) = £ {^MqMsMj^IMt)^ ® |lM s ) Slight } 

M S M Q 

x ^ Q (r 1; r 2 ,r 12 )®|Q;M Q ) (17) 
\X* Q ;\J=^,Mj;M T ) = £ ( ^l||MgM s Mj){|lM r )i ® |lMs)s |Ught } 

Ms Mq 

x ^(n.ra,^)® |<5;M Q ) (18) 

with Mr and Mq, the baryon isospin and heavy quark spin third components, respectively. 

As in the Aq case, and based on the HQS predictions, here we have also neglected components constructed out 
of the spin-singlet (|Q0)<; U , t ) and anti-symmetric spatial wave functions, since they are suppressed by powers of 
1/toq. HQS leads to similar simplifications for the rest of baryons studied below and we will omit any further 
comment in what follows. 

S— type baryons: I = h , Slight = 

|H Q ;J=i Mj;M T ) = ^={\ls)^f s Q ( ri ,r 2 ,r 12 ) - \sl)*l Q ( ri ,r 2 ,r 12 )} ® |00> Sught ® \Q; Mj) (19) 

where the isospin third component of the baryon, Mr, is that of the light quark / (1/2 or —1/2 for the u or the 
d quark, respectively). 

H' and S*—type baryons: I — j. Slight = 1- 

|5 Q ;J=i Mj;M T ) = J2 {\^\\M Q M s Mj)±= 

m s m q ^ 

X {|^)< s b (n,r 2 ,r 12 ) + \sl)*f (n,ra,r 13 )} ® |lM g ) Slight ® \Q;M Q ) (20) 
\E* Q ; J = 1 Mj; M T ) = £ (lil\M Q M a Mj)±= 

M S M Q VZ 

x {M*?(n, ra.ru) + |sO*?(n,r 2 ,n 2 )} ® |lM s ) Slight ® |Q;M Q ) (21) 



where the isospin third component of the baryon, Mt, is that of the light quark 
• and Q* — type baryons: 1 = 0, Slight = 1- 

P Q ;J = \,Mj ) = E dl^lMoMsMjJIlM^^^Cri.ra.ru) \Q;M Q ) (22) 

M s Mq 

|n^J=|,Mj) - Yj (^l\M Q M s Mj)\lM s ) Snsht ^ Q (n,r 2 , ru )\®\Q-,M Q ) (23) 

m s m q 
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where 4^ S Q (n, r 2 , r 12 ) = * 5S Q (r 2 , n, r X2 ). 

The spatial wave function, 4"^? , will be determined by the variational principle: <5 | iJ lnt | J3q } =0. For simplicity, 
we will use a family of functions with free parameters, which will be determined by minimizing the mass of the baryon. 
The OGE has a regularized delta function, which provides strong short range repulsions/attractions, specially between 
the light quarks 4 . A similar situation appears in the context of double A hypernuclei |3 1|. which in a very good 
approximation can be also reduced to a three body problem. In this context, it was shown that the use of a Jastrow- 
type functional form for the spatial wave functions leads to excellent results, and to notably simpler wave functions 
than other functional forms, as for instance a series of standard Hylleraas type wave functions |35|. Thus, we take 



Ori,r 2 ,r 12 ) = NF B ^(r 12 )^(n)^,(r 2 ) 



(24) 



where N is a constant, which is determined from normalization 5 



ri 



r 2 



* OT ?(ri,r 2 ,r 12 ) 



Sir* 



+ OC 



+ 30 



dr 2 



+i 



(25) 



where [i is the cosine of the angle between the vectors ri and r 2 , being ri 2 = ( r\ + r\ — 2r\r 2 [i) 2 . For simplicity, we 
do not entirely determine the functions (j>® and <j>®, from the variational principle, but we rather fix the bulk of these 

functions to the s— wave ground states (tfi{ =q q i) of the single particle hamiltonians, h^L qql , defined in Eq. (JSJ), and 
modify their behavior at large distances. Thus, we take 



<^(n) = (i + a g n)^(n) 

€(r 2 ) = (l+tvra)¥>S(ra) 



(26) 



with only one (two) free parameter for a 11 or ss (Is) baryon light quark content. Besides, we construct the light-light 
correlation function, F Bq , from a linear combination of gaussians, 



-«(ria+<if)' 



, ai = 1 



f B Hr 12 ) 



!_ e ~cr 12 if Vq B (r 12 = 0) >0 
1 (c->+oc) if V B ,(r 12 =0)<0 



(27) 
(28) 



where V B q , denotes the light-light interaction projected onto the spin and isospin quantum numbers of the baryon B. 
The correlation function, F Bq , should vanish at large distances because of the confinement potential 6 . The function 
J Bq vanishes at the origin in those cases where V q B , is highly repulsive a short distances (this is the case for instance for 
the ALlx interaction and Sii g ht = 1, i-e. S,S*— type baryons). Note that, one of the aj parameters can be absorbed 
into the normalization constant N, thus besides the c and ct q%q i parameters, there are a total of eleven (we fix a\ = 1) 
free parameters to be determined by the variational principle and the mass of the baryon is just the expected value 
of the intrinsic hamiltonian. Note that all interactions of Subsect. Ill Bl are diagonal in the flavor space, except for 
the pion exchange term of the AL1% one, and that the following relationships are of interest to compute the kinetic 
energy expected values 



Vi 

v 2 



4>i{ri)4> 2 {r 2 )F(ri 2 ) 
0i (ri)(/> 2 (r 2 )F(ri 2 ) 



fxd4>i{r x ) r x -r 2 ,dF(r 12 ) 
-4> 2 (r 2 )F(r 12 ) H 4>i(ri)(p 2 (r 2 ) 



ri dr\ 



ry 2 



dr\ 2 
dF(r 12 ) 



r 2 , , ,# 2 (r 2 ) r 1 ~r 2 

— (Pi(ri) — 1 F(n 2 ) (Pi(ri)(p 2 (r 2 ) 

r 2 dr 2 r 12 dr 12 



(29) 



4 The size of this term of the OGE interaction is proportional to the inverse of the masses of the involved quarks. 

5 Note that, the Jastrow form assumes a factorization of the wave function in three terms, each of them depends only on one of the 
involved variables ri, r2 and r\2- Dependences of the type ri — T2, mentioned in the discussion of Eq. 1161 . cannot be accommodated 
with factorizable functions of this type. 

6 The confinement potential is also responsible for the non— vanishing values of the parameters a q and a' q in Eq. 1261 . Indeed, if at large 

distances the light— light interaction vanishes (V^,(ri2 ^ 1) — » 0) then the product of single particle wave functions ip^(ri) </^, (r'2) will 
provide the correct spatial dependence, at large distances, of the solution of the three body problem. 
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The expected value of the Hughes-Eckart contribution (Vi • V2) to the kinetic energy is readily obtained from the 
above equations after integrating by parts. 

Given the spatial part of the baryon wave function, ^'^ < ? (^i, f 2 , r 12 ), the probability, Vi, of finding each of the two 
light quarks with angular momentum I, referred to the heavy quark, and coupled to L = is given by 



r + 00 p+OO p+l 

Vi = 4tt 2 (2Z + 1) / dnrf / dr 2 r 2 2 / dfj, Pl(ji)**5 (n, r 2 , r X2 ) 

Jo Jo J-i 

where Pi is the Legendre Polynomial of order I . 



(30) 



IV. STATIC PROPERTIES: MASS AND CHARGE DENSITIES AND ORBITAL MAGNETIC 

MOMENTS 



Since the charge operator is diagonal in the spin-flavor space, the baryon charge density at the point P (coordinate 
vector f in the CM frame, see Fig. ^) is given by: 



d^R d 3 nd 3 r 2 



p iPauR ,, 

' ' He 



* g< ?(ri,ra,ria) {e Q S 3 (r - y h ) + e q 6 3 (r- y t ) + e q ,5 3 {r- y 2 )} 



= J d 3 r 1 d 3 r 2 



*f?(ri,ra,ria) " {e Q S 3 {r - y h ) + e q S 3 {r~ y x ) + e q ,8 3 {r~ y 2 )} 



Pf Q (r)\ Q +p^(r)\ q + p^(r)\ ql 



(31) 



where <$(• • •) and Pcm are the three-dimensional Dirac's delta and the total baryon momentum respectively, V is 
the interacting volume 7 , e q . q \Q are the quark charges, in proton charge units (e), and from Fig. ^ we have 8 , yn = 
— [m q f\ + m q ir 2 ) /M, yi = yk+ri and y 2 = yh+r 2 . The charge density is spherically symmetric (pe Q (f) = Pe Q (|f |)), 
since the wave function only depend on scalars (r 1; r 2 and fx -f 2 ) for L = baryons. The charge form factor is defined 
as usual 

(32) 



F? Q m= J d 3 re^p^(r) 
and it only depends on \q |. The charge mean square radii are defined 

p p + oc 

( r2 )e Q = J d 3 rr 2 pfi(r) = 4tt J drr 4 p^{r) 



(33) 



The baryon mass density, pm ' (r) is readily obtained from Eq. (|31|1 with the obvious substitutions (e q ,e q / ,cq) 
(m q /M,m q ,/M,m Q /M). 



where we have normalized pm Q (r) to 1. The mass mean square radii are defined 



(r 2 )^ = I d 3 rr 2 p^(r) = 4tt J drr 4 p%> (r) 



(34) 



(35) 



The intrinsic orbital magnetic moment is defined in terms of the velocities Vh.,1,2 of the quarks Q, q and q 1 , with respect 
to the position of the CM, and it reads 



d*R d 3 nd 3 r 2 



-iPcmR 



(C< 



ri,r 2 ,r 12 . 



[ e Q 



I 2mQ 



(Vh x m Q v yh ) 



G & 1 I 6* 



-*f 9 ?(ri,r 2 ,r 12 ) 



(36) 



It cancels out with J d 3 R. 

There exists the obvious restriction mQy^ + m q y\ + m q /y2 = 0. 
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with 9 m q v yi = — iVi, m q iv y2 = — iV 2 and mQV yh = i (\7i + V 2 ^ ■ Since the wave function only depends on scalars 

for a L — baryon, the intrinsic orbital magnetic moment vanishes. Note that the intrinsic orbital magnetic moment 
is zero despite the quark pairs (Q, q) and (Q, <?') are not in relative s— waves, as mentioned above in the discussion of 
Eq. I|30|l . Thus, the magnetic moment of the baryon is entirely due to the spin contribution and since we are adopting 
HQS wave functions where the light degrees of freedom spin wave function factorizes, the magnetic moment should 
coincide to that obtained in the naive quark-model treatment of the baryon. 

Finally and as a further test of our variational wave-functions, we compute some quantities which in Ref. |28l | are 
denoted by Pi(0) and called wave function at the origin for the (jk) pair, where (i,j,k) is a cyclic permutation of 
(Q, 9i ?')• I n this reference the following Jacobi coordinates 



SQ(fl,f 2 ) 

s q {ri,r 2 ) 



n - r 2 
<V* 2 , 



s q '{r 1.^2) = a q >n, a q > 



l m Q (m q + m ql ) 
m q {mQ + m q ') 

/ m Q {m q + m q ,) 
y m q '(iriQ + m q ) 



(37) 



and correlation functions 
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vT Q (s) 



d 3 r 1 d 3 r 2 



#J?(n,r2,ri 2 ) S(s- |sf(fi,f 2 )|) , i = Q,q,q' 



(38) 



are introduced. Thus, the wave functions at the origin are defined by [28 

r +00 



0) 



B «(s = 0) 



V* Q (s = 0) 



16tt 2 
16tt 2 

16tt 2 



dtr 



+ 00 



dW 



-HCxD 



dtr 



99' 1 



t,r 2 = t,ru = 
= t,r 2 = 0,n 2 
= 0,r 2 = i,ri 2 



0) 



i) 



t) 



(39) 



V. RESULTS 



In Table 11111 we present variational results for the charmed and bottom baryon masses obtained with each of the 
six quark-quark interactions considered in this work. We use an integration step of 10 -2 fm and we estimate in 
0.5 MeV the numerical error of the variational masses given in the table. On the other hand, by including a higher 
number of gaussians 11 in the light-light correlation function (Eq. Q27p). or quadratic or cubic terms in the one body 
wave-functions (Eq. (|26|l ). the variational energies cannot be lowered in more than 0.5 MeV. The corresponding wave 
functions can be reconstructed from the parameters given in Tables Ixl and fxTI in the Appendix. We use a Simplex 
algorithm to find out the position of the minima [36j and typically for each baryon the minimization procedure, 
including the construction of the wave function, takes around twenty minutes in a 2 GHz Pentium IV processor. 

Besides, in Table UTTl we compare, when possible, with the Faddeev results of Ref. [28|. We find an excellent 
agreement between our variational results and those of Ref. [2^|. In most cases we find discrepancies smaller 12 than 3 
MeV, and only in three cases (BD il c , and AL2 Af, baryons) Faddeev and variational masses differ in more than 5 
MeV. Faddeev masses are usually smaller than variational ones, but in some cases the variational approach presented 
here provides lower values, and thus better estimates of the baryon masses (see f.i. AL2 Aj, where we find a mass of 



9 Note that the classical kinetic energy has a term on v yi ■ v y2 and then the operator m q v yi is not proportional to —i V W1 , but it is rather 
given by m q v yi = (— i V yi ) - V y2 ) = (-i Vi). Similarly m q ,v y2 = (-i V2). 

These correlation functions satisfy J dss 2 r] i Q (s) = 1. 

1 In the Appendix (Table ITxl . we illustrate, for the A c baryon with an AL1 inter-quark interaction, the dependence of the VAR results 
on the number of gaussians. 

2 Such differences are certainly small fractions of the baryon binding energies. 
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Charm Bottom 



Baryon 


BD 


AL2 API AP2 AL1 AL1 X 


Exp 


BD 


AL2 API AP2 AL1 AL1* 


Exp 




FAD 


2321 


2302 


2311 


2313 


2296 


— 




5660 


5637 5659 


5643 


5643 


— 




A 


VAR 


ooon 
ZoZ[) 


2305 


2308 


2310 


2295 


2140 


OOOC _l_ 1 


ODDO 


5629 


5659 


5643 


5643 


5475 


obz4 ± y 




FAD 


OA Q/1 
Z4y4 


2474 


2484 


2491 


2466 






OoDO 


5844 


5868 


5859 


5849 






£ 


VAR 


9/1 Qft 
z4yo 


2479 


2484 


2491 


2469 


2457 


0/1 co -1- 1 
Z4DZ zt 1 


Do ( 1 


5844 


5870 


5861 


5851 


5841 


^77nt -1- 70 
ii U ± (U 


E* 


VAR 


2570 


2559 


2560 


2571 


2548 


2535 


2518 =b 2 


5900 


5874 


5899 


5893 


5882 


5869 


5780' ± 70 




FAD 


2502 


2480 


2485 


2488 


2473 






croon 


5803 


5819 


5805 


5808 








VAR 


2506 


2480 


2485 


2488 


2474 




OARQ -4- Q 


croon 
OooZ 


5800 


5820 


5807 


5808 




c'varxt _|_ 7n 
/ OU ± / u 




VAR 


2606 


2584 


2587 2592 


2578 




2576 ± 3 


5963 


5939 


5959 


5946 


5946 




5900!^ ± 70 




VAR 


2671 


2661 


2664 


2672 


2655 




2646 ± 2 


5989 


5970 


5987 


5979 


5975 




5900!^ ± 80 




FAD 


2707 2686 


2678 


2682 


2678 






6048 


6029 


6033 


6020 


6035 








VAR 


2713 


2683 


2682 


2684 


2681 




2698 ± 3 


6050 


6030 


6036 


6023 


6033 




5990 1 " ± 70 


n* 


VAR 


2770 


2759 


2759 


2765 


2755 




2660 1 " ± 80 


6074 


6061 


6068 


6056 


6063 




60oo T ± 70 



TABLE III: Masses (in MeV) of the ground states of different baryons with the six potentials considered in this work. Variational 
and Faddeev results are given in the rows denoted by VAR and FAD, respectively. The latter ones are taken from Ref. |28j, but 
the effect of the three body forces considered in that reference, V123 = constant /m q m q imQ, has been eliminated. Experimental 
masses are taken from Ref. [To| when possible, and the lattice QCD estimates of Ref. [l2T| are quoted for those cases (f) where 
the experimental masses are not known. 



around 8 MeV smaller than that quoted in Ref. [28]] )• Having in mind that the Faddeev masses suffer from numerical 
uncertainties of the order of 5 MeV [2^] , we can safely conclude that our variational approach reproduces the ground 
state masses quoted in Ref. Thus, corrections to the leading HQS wave functions (Eq. (fTo")) and Eqs. l(T7l) - (|25|) ') 
should provide quite small (let us say, smaller than 5 MeV) changes in ground state masses, even for the charm sector. 

We have also computed the E*, H', S* and fl* ground state masses, which were not evaluated in the work of 
Ref. [2j|. Note that, the a q ■ Gqjm q mQ OGE term breaks 13 the mass degeneration between the Slight triplet states 
(£-£*, S'-S* and O-O*). The five phenomenological interactions studied in Ref. [2^ lead to ground state masses 
which compare reasonably well with the existing experimental data and lattice QCD mass estimates for all baryons. 

Besides, we see that while the AL1 and ALlx interactions lead to similar masses (differences of the order of 10 MeV) 
for the Eq and Eg baryons, the chirally inspired potential, ALlx, predicts masses for the Ag baryons of about 155 
(charm sector) and 168 (bottom sector) MeV smaller than those obtained from the phenomenological AL1 interaction. 
In terms of binding energies (Bq = M [Aq] — mg — m u — rrid), we have relative changes as large as 90%, 



B AL1 _ oALlx B AL1 _ rAL1 x 



B ALl 



0.9, 



0.7 (40) 



Such big changes are due to the different behaviour of the One Pion Exchange (OPE) potential in the light quark- 
quark and quark-antiquark sectors. In the first case we have, from Eq. i|12[l . the operator (diUj) (fifj) which provides 
a factor 9 (1) for the light quantum numbers of the Aq (Eq or Eg) baryon. The spin-isospin operator turns out to be 
— {(JitJj) (TiTj) for the case of quark-antiquark interactions, and thus one gets a factor 3 (—1) for the light quantum 
numbers of the tt (p) meson. In averaged, both potentials, AL1 and AL1%, are strongly attractive (moderately 
repulsive) for the 7r (p) quantum numbers, and have been adjusted to provide the same ground state pion (rho) mass. 
For baryons, the AL1 potential gets reduced its strength by a factor of two (because of the V^ q = V qq /2 prescription, 

coming from a XiXj color dependence), however this is not the case for the light sector of the AL1\ potential. In this 
latter case, the OGE and confinement terms follow the former prescription, but while the scalar exchange remains 
unaltered, the OPE term is three times more attractive for the Aq baryon than in the pion. The resulting effect 
is that the light sector of the ALlx interaction in the Aq baryon gets reduced its strength by a factor significantly 
smaller than 2. The OPE exchange plays a minor role in the spin triplet states 14 (p, Eq and Eq), and that explains 
why the AL1 and AL1% potentials predict similar Eq and Eg masses. 



For s-wave baryons, the spin operator a q ■ <tq vanishes for Slight = wave functions and for spin triplet states, Slight = 1 5 takes the 
values 1 and —2, for baryon total angular momentum J = § and i, respectively. 

For instance, while the OPE is responsible for around 150 MeV of the binding energy of the pion, it contributes only about 8 MeV to 
the rho meson mass. 
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Rat'vnn 


Charm 

RD AT.9 API AP9 AT.1 AT,1v 


Bottom 

RD AP9 API AP9 API AT.1v 


FAD 
A VAR 


0.097 0.101 0.109 0.103 0.104 — 
0.100 0.101 0.108 0.106 0.106 0.101 


0.043 0.044 0.046 0.044 0.045 — 
0.045 0.043 0.048 0.045 0.045 0.042 


FAD 

E VAR 


0.111 0.117 0.128 0.121 0.121 — 

Alio n no n i 07 Alio n i oo n i oo 

U.iiz U.iio U.lzY U.iio U.L26 U.izo 


0.051 0.053 0.056 0.053 0.054 — 

n ncri n no n nnc n no n ak7 n nt^n 

0.051 0.053 0.055 0.053 0.057 0.059 


E VAR 


0.122 0.131 0.138 0.131 0.135 0.139 


r\ ata r\ at/ 1 n ato n r\t~ f r\ f\rr\ r\ aai 

0.052 0.056 0.058 0.056 0.060 0.061 


FAD 
S VAR 


0.097 0.100 0.107 0.101 0.104 - 

n nnn n i no nun n i n^i n i nt^ 

0.099 0.103 0.110 0.104 0.105 — 


0.045 0.046 0.048 0.046 0.048 - 

n n at r\ r\ a a nncn r\ r\ a*~7 n n ac\ 

0.047 0.046 0.050 0.047 0.049 — 


WAT} 


n i (it; n nn n i oo n n 1 n 1 1 n 
0.106 U.110 U.122 U.lll U.119 — 


O.Uoo O.Uoo 0.U59 0.U52 0.U6U — 


H* VAR 


0.115 0.123 0.133 0.126 0.123 - 


0.052 0.057 0.060 0.054 0.059 - 


FAD 
n VAR 


0.100 0.103 0.110 0.102 0.108 - 
0.096 0.101 0.110 0.103 0.108 - 


0.050 0.052 0.055 0.051 0.054 - 
0.050 0.054 0.057 0.053 0.057 - 


fi* VAR 


0.107 0.117 0.122 0.113 0.120 - 


0.051 0.055 0.061 0.056 0.059 - 



TABLE IV: Mass mean square radii in mi 2 units, see Eq. Ij35^ . Variational and Faddeev |2Sj results are denoted by VAR and 
FAD, respectively. 



_L> oil y Uil 


Charm 

RD AT,° API AP9 AT-1 AT,1v 


Bottom 

RD AT, 9 API AP9 API AT,1v 

Uly jVJJij ill _1_ ill £i AlJl Ix. 1— J _L X 


FAD 

A+, A° b VAR 


0.117 0.124 0.147 0.139 0.129 

0.120 0.125 0.145 0.142 0.131 0.124 


0.115 0.123 0.148 0.140 0.128 

0.121 0.122 0.152 0.143 0.127 0.119 


FAD 

E°, E" VAR 


-0.224 -0.245 -0.304 -0.287 -0.256 
-0.226 -0.248 -0.299 -0.279 -0.261 -0.260 


-0.280 -0.305 -0.369 -0.347 -0.318 
-0.279 -0.304 -0.357 -0.341 -0.332 -0.345 


FAD 

E+, E° VAR 


0.134 0.145 0.174 0.164 0.151 

0.135 0.147 0.172 0.160 0.153 0.154 


0.138 0.151 0.183 0.172 0.157 

0.138 0.150 0.176 0.169 0.164 0.171 


FAD 

E++, E+ VAR 


0.494 0.535 0.652 0.615 0.557 

0.497 0.541 0.643 0.599 0.568 0.567 


0.555 0.607 0.734 0.692 0.633 

0.555 0.603 0.710 0.679 0.659 0.687 


E* u , E*~ VAR 


-0.246 -0.274 -0.323 -0.307 -0.283 -0.293 


-0.289 -0.322 -0.376 -0.364 -0.349 -0.358 


Er, Sr VAR 


0.148 0.163 0.187 0.177 0.168 0.173 


0.143 0.159 0.186 0.180 0.172 0.177 


E* ++ , K + VAR 


0.541 0.599 0.697 0.661 0.619 0.639 


0.574 0.639 0.747 0.725 0.694 0.711 


FAD 

H°, H" VAR 


-0.145 -0.154 -0.183 -0.171 -0.161 
-0.148 -0.156 -0.182 -0.174 -0.163 


-0.193 -0.203 -0.234 -0.219 -0.212 
-0.198 -0.201 -0.240 -0.220 -0.213 


FAD 

S+, H° VAR 


0.160 0.170 0.206 0.196 0.177 
0.158 0.170 0.198 0.196 0.179 


0.151 0.161 0.197 0.187 0.168 
0.150 0.151 0.189 0.180 0.158 


~™, E'f VAR 


-0.164 -0.179 -0.213 -0.196 -0.192 


-0.228 -0.235 -0.285 -0.248 -0.267 


H^, VAR 


0.174 0.199 0.234 0.218 0.207 


0.184 0.188 0.234 0.209 0.214 


H^ u , H*" VAR 


-0.177 -0.196 -0.232 -0.217 -0.198 


-0.224 -0.250 -0.287 -0.259 -0.266 


H c * + , H b * u VAR 


0.190 0.218 0.254 0.241 0.213 


0.179 0.201 0.236 0.221 0.212 


FAD 

Q°, 0,7 VAR 


-0.111 -0.117 -0.131 -0.120 -0.124 
-0.108 -0.115 -0.131 -0.120 -0.124 


-0.164 -0.173 -0.191 -0.174 -0.183 
-0.163 -0.177 -0.196 -0.180 -0.189 


n* u , n*- VAR 


-0.118 -0.132 -0.143 -0.131 -0.138 


-0.167 -0.182 -0.210 -0.189 -0.196 



TABLE V: Charge mean square radii in efm 2 units, see Eq. I|33|l . Variational and Faddeev |2Sj results are denoted by VAR 
and FAD, respectively. 

We have also computed mass ( Table HVj) and charge (Table W} mean square radii, charge form factors (Figs.HH3) 
and quark mass densities (Figs. 0H3)i f° r au baryons and the six inter-quark interactions considered in this work. 
In Tables IIVI and when possible, we compare our variational results with those quoted in Ref. j2^| and find an 
excellent agreement, with tiny differences which at most are of the order of 3-4%. Differences between AL1 and ALlx 
potentials predictions, in Tables llVl and IVl and in Figs. [2HS1 are a ls° quite small and it will be quite difficult to use 
those to experimentally disentangle among both interactions. 

To further test the goodness of our variational baryon wave functions, we have followed Ref. [2^ and we have 
also computed the so called wave functions at the origin for the different quark pairs inside the heavy baryon (see 
Eq. I|39I) L It seems [2^ that the absolute value of this quantity depends dramatically on the numerical procedure 
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used to solve the three body problem 15 and then is a perfect observable to check the validity of our approach. Results 
for the different inter-quark interactions and heavy baryons can be found in Table IV1I where, when possible, we 
compare our variational results with those quoted in Ref. pflj . as well. For A—, S— and Q— type baryons, we find 
that our variational results and those obtained from the Faddeev calculation of Ref. [28j, agree 16 within a 10 or 20%. 
This constitutes a remarkable result, since, as it is mentioned in Ref. |2q. it is not surprising to have approximate 
treatments of the three body problem giving values for the wave function at the origin off by almost an order of 
magnitude. 

However, variational and Faddeev results are in total disagreement for the H baryons (discrepancies of about a 
factor three or larger ). W e have no explanation for this fact, beyond that saying that we might have misunderstood 
the definition of Ref. [28| of the wave functions at the origin and thus Eq. I|39ll is wrong or that our S— wave functions 
at the origin are wrong. We very much doubt this latter possibility, and we would like to point out that the Faddeev 
values for the S— baryons are abnormally big (about a factor four), when compared to those quoted for the rest of 
the baryons A, £ and Q. We are not aware of any reason to support/explain this fact and indeed we do not observe 
such a situation in the variational results. 

To finish this section we would like to devote a few words to the multipolar content of the variational wave functions. 
In Tables I VIII and IVIIII we show the probability Vi , defined in Eq. (|30|l , of finding each of the two light-quarks with 
angular momentum I and coupled to L = for the A and S baryons. Results obtained with two different interactions 
(AL1 and AL1%) are compiled, but the gross features of the multipolar decomposition do not depend significantly on 
the selected interaction. Only multipoles up to the waves I = 3 or I = 4 significantly contribute 17 , being I — the 
dominant one. 



VI. CONCLUSIONS AND OUTLOOK 



We have developed a variational scheme to describe the charm and bottom baryons compiled in Table [I] We 
compute baryon masses and wave functions from the variational principle applied to a family of Jastrow type functions 
constrained by HQS. For several inter-quark interactions and baryons, we reproduce the Faddeev results of Ref. j2^ 
for masses, charge and mass radii, and wave functions at the origin. Thus, we conclude that the finite heavy quark 
mass corrections to the wave-functions should be small, even for charmed baryons. Besides of giving results for 
baryons not studied in Ref. [22j (£*, £', S* and ft*), we provide, in Tables IXl and IXII in the Appendix, wave functions, 
parameterized in a simple manner, for all baryons compiled in Table [I] and six different inter-quark interactions. 
Thanks to HQS, the baryon wave functions are significantly simpler and more manageable than those obtained from 
Faddeev calculations and they can be used by the community to compute further observables, without having to solve 
the involved three body problem. Hence, we have shown in this context, once more, that HQS is a useful tool to 
understand the bottom and charm physics. We have also worked out the predictions for the strangeness-less baryons 
of the SU(2) chirally inspired quark-quark interactions developed in Ref. j2^,|2!|. The chirally inspired potential, 
ALlx, predicts masses for the Aq baryons of about 155 (charm sector) and 168 (bottom sector) MeV smaller than 
those obtained from the phenomenological AL1 interaction. Assuming that the inclusion of a pattern of SCSB should 
lead to a more theoretically founded scenario to address hadron spectroscopy, the findings of this work might hint 
the existence of sizeable three quarks forces with the Aq quantum numbers. New SU(3) chirally inspired quark-quark 
interactions, valid also for the s quark (exchanges of eta and kaon mesons are also introduced) are being developed 
by the Salamanca group |37|. It would be interesting to work out the predictions for the strangeness baryons with 
this new set of chirally inspired interactions, once the potential becomes available. 

Using the semi-analytical wave functions found in this work, we plan to compute the scmilcptonic decay of the 
bottom baryons into the charmed ones. The mixing parameter \V c b\ of the Cabibbo-Kobayashi-Maskawa matrix can 
be determined from the study of these processes. Furthermore, the Isgur-Wise form-factor, which is a universal 



Indeed, the wave function at the origin is especially sensitive to short-range correlations which are not precisely taken into account by 
approximate treatments. 

Note that for S baryons and when the BD interaction is used, there exists a serious discrepancy between Faddeev and variational 
predictions for the 11 pair wave function at the origin. Indeed, the Faddeev prediction is non-zero, while it is zero within the variational 
approach. This is because our variational scheme assumes a pure Slight = 1 configuration for this baryon. It turns out, that the spin 
triplet light quark-light quark BD potential is infinitely repulsive for r±2 = and conversely the variational wav e function vanishes at 
the origin (because of the correlation function / s< 3 , in Eq. 1281 1. However in the Faddeev calculation of Ref. l2St . the HQS suppressed 
•Slight = component of the E wave function is kept, and since for the spin singlet configuration the BD potential becomes attractive, 
a non- vanishing wave function at the origin for the 11 pair is found in Ref. l2Sl . 

Note that we show probabilities: the components of the wave function for each value of I are given, up to a sign, by the squared root of 
the numbers presented in the tables. 
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Baryon Pair 



Charm 
BD AL2 API AP2 AL1 

15.0 12.1 

18.5 14.0 

7.7 6.1 

7.0 5.6 



AL1 X 



BD AL2 



Bottom 
API AP2 



AL1 AL1 X 



A 



11 
11 

Ql 
Ql 



FAD 
VAR 
FAD 
VAR 



23.1 
22.8 
8.8 
7.6 



11.9 14.8 

13.9 18.6 

6.4 7.4 

5.5 6.5 



54.0 
6.0 



23.3 
19.5 
8.8 
7.5 



15.2 
19.1 
7.7 
7.9 



12.2 
14.4 
6.1 
4.9 



12.1 
14.8 
6.4 
5.6 



14.9 
18.3 
7.4 
6.9 



74.7 



7.0 



11 
11 
Ql 
Ql 



FAD 
VAR 
FAD 
VAR 



6.9 
0.0 
8.2 
7.3 



6.6 5.0 
6.8 5.2 

6.7 5.3 
6.0 4.6 



5.6 
6.2 
5.5 
5.2 



6.1 
5.7 

6.5 
5.9 



0.0 



5.5 



6.8 
0.0 
7.1 
6.7 



6.4 
6.9 
5.9 
6.1 



4.9 
5.5 
4.7 
3.9 



5.4 
5.3 
4.9 
4.5 



6.0 
6.7 
5.7 
4.8 



0.0 



4.2 



E* 



11 

Ql 



VAR 
VAR 



0.0 
5.7 



6.6 
4.3 



4.7 
3.5 



5.3 
3.8 



5.9 
4.4 



0.0 
4.1 



0.0 
6.4 



6.3 
5.2 



5.1 

3.3 



5.2 
3.5 



6.4 
4.0 



0.0 

3.9 



Is 
Is 

Ql 
Ql 
Qs 
Qs 



FAD 
VAR 
FAD 
VAR 
FAD 
VAR 



70.5 
21.1 
40.0 
11.7 
30.8 
8.8 



54.6 
21.9 
37.0 
9.5 
27.7 
7.5 



51.0 
17.0 
34.4 
8.3 
25.3 
6.1 



52.1 
18.1 
37.0 
8.7 
26.4 
6.8 



53.5 
19.3 
35.5 
9.1 
27.1 
7.8 



93.1 
23.8 
51.4 
11.6 
39.9 
8.7 



72.0 
20.4 
48.1 
11.9 
35.5 
8.3 



66.8 


68.4 


70.2 


15.6 


17.2 


20.2 


44.6 


48.0 


46.1 


8.7 


10.2 


10.0 


32.4 


33.4 


34.9 


6.1 


7.2 


7.3 



Is 
Ql 
Qs 



VAR 
VAR 
VAR 



12.5 
11.3 
8.7 



10.3 8.4 
9.5 7.1 
8.2 5.8 



9.2 9.6 
9.0 8.2 
6.9 6.5 



11.9 
9.8 
8.2 



10.1 
9.5 
7.3 



8.2 
5.8 
4.7 



9.1 
8.6 
6.3 



9.2 
6.7 
5.5 



Is 

Ql 
Qs 



VAR 
VAR 
VAR 



12.1 
8.5 
6.7 



8.8 6.8 
6.4 5.8 
5.7 4.8 



7.9 8.9 
6.0 7.7 
4.6 6.1 



12.1 
10.4 
8.4 



9.3 
8.0 
6.3 



7.9 
5.7 
4.7 



8.4 
7.7 
6.0 



9.4 
6.8 
5.5 



ss 
ss 

Qs 
Qs 



FAD 
VAR 
FAD 
VAR 



16.5 
17.7 
17.8 
17.0 



15.5 13.8 

17.2 14.8 

15.4 14.4 

14.5 12.3 



15.3 14.6 

16.9 15.7 

15.0 15.1 

13.5 13.0 



16.7 
19.9 
16.0 
15.3 



15.6 
15.8 
13.8 
12.5 



13.7 
13.1 
12.8 
10.8 



15.3 
16.9 
13.3 
11.6 



14.6 
15.5 
13.5 
11.1 



ss 
Qs 



VAR 
VAR 



18.7 
12.6 



13.5 
9.9 



13.5 
9.4 



14.9 
10.2 



13.0 
9.7 



18.4 
15.0 



15.7 
11.7 



12.5 
9.4 



15.0 
10.7 



14.3 
10.4 



TABLE VI: Wave function at the origin, rjf Q (0) in fm 3 units, see Eq. 139^ . In the second column, 11 and ss pairs stand for 



rig® (0) and r/Q Q (0). Besides for strange [strangeness-less ] baryons, the pairs Ql and Qs stand for 7^,5^(0) 



and 7/ 9=; (0). Variational and Faddeev |28j| results are denoted by VAR and FAD, respectively. For the E and E* baryons, our 

b" 

variational scheme assumes a pure Slight = 1 configuration. The VAR estimates for t]q Q (0), and BD and ALlx interactions, 
vanish for these baryons, since the spin triplet light quark-light quark forces are infinitely repulsive at r\2 = 0, for these 
potentials. 



1 


Charm 

1 2 3 4 5 6 


Bottom 

1 2 3 4 5 6 


ALlx-Vi 


0.75495 0.18706 0.04066 0.01103 0.00368 0.00143 0.00062 
0.7550 0.9420 0.9827 0.9937 0.9974 0.9988 0.9994 


0.79079 0.16186 0.03237 0.00894 0.00324 0.00139 0.00067 
0.7908 0.9526 0.9850 0.9940 0.9972 0.9986 0.9993 


ALl-Vi 
J2k=n 


0.87389 0.10972 0.01346 0.00227 0.00052 0.00015 0.00005 
0.8739 0.9836 0.9971 0.9993 0.9999 1.0000 1.0000 


0.87368 0.11207 0.01221 0.00169 0.00032 0.00008 0.00003 
0.8737 0.9857 0.9980 0.9996 1.0000 1.0000 1.0000 



TABLE VII: Probabilities Vi (defined in Eq. iJHSJ) for several waves. Results were obtained for the A baryon and AL1 and 
ALlx inter-quark interaction. The errors are always less than one unit in the last digit. 



/ 


Charm 

1 2 3 4 5 6 


Bottom 

1 2 3 4 5 6 


AL1 X -Vi 


0.95074 0.04711 0.00210 0.00014 0.00002 0.00000 0.00000 
0.9507 0.9978 0.9999 1.0000 1.0000 1.0000 1.0000 


0.95010 0.04777 0.00204 0.00013 0.00002 0.00000 0.00000 
0.9501 0.9979 0.9999 1.0000 1.0000 1.0000 1.0000 


KLl-Vi 


0.96542 0.03377 0.00091 0.00003 0.00000 0.00000 0.00000 
0.9654 0.9992 1.0000 1.0000 1.0000 1.0000 1.0000 


0.96109 0.03734 0.00153 0.00011 0.00001 0.00000 0.00000 
0.9611 0.9984 1.0000 1.0000 1.0000 1.0000 1.0000 



TABLE VIII: Same as Table IVTTI for the E baryon. 
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FIG. 2: Charge form factors (Eq. 1321 ) of strangeness-less charmed and bottom baryons obtained with AL1, ALlx interactions. The only 
appreciable difference, as it is also the case for the spectrum, appears for the A c and A& baryons. Note that, ^(0) determines the charge 
of the baryon. 




FIG. 3: Charge form factors (Eq. 1321 ) of strangeness charmed and bottom baryons obtained with AL1 and BD quark-quark interactions. 
The charge of the baryon is given, in units of the proton charge, by the value of the form factor at the origin. 
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FIG. 4: Quark mass densities [fm 3 1 , p„^ (r) I andp„ <3 (r*)| (Eg. 1341 V of strangeness-less charmed and bottom barvons obtained 

\Q — C,b \q — l 

with ALl and ALlx quark— quark interactions. 
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IN UlllUttI Ul ^cLLlholdllo 


A Ma«s rMeVl 
iV c IVicLab llvlc V 1 


W [tm j 


n Ac ('n^ rfm _3 i 
V c i (UJ [tm J 


I 


2302.1 


16.4 




2 


2295.0 


17.0 


5.7 


3 


2294.8 


18.6 


6.2 


4 


2294.6 


18.6 


6.5 


FAD 


2296 


14.8 


7.4 



TABLE IX: A c baryon mass and wave functions at the origin (see Eq. £2)) as a function of the number of gaussians included 
in the light-light correlation function. Results have been obtained with the AL1 inter-quark interaction and for comparison, 
the Faddeev results of Ref.|2g] (denoted by FAD) have been also compiled. 



function that governs all the exclusive b — ► c semileptonic decays in the tuq — > oo limit, could be obtained and more 
importantly the size of the \jrriQ corrections in these decays would be also estimated. We will compare our results 
to those obtained in lattice QCD simulations, both from the meson and baryon sectors |3Sj . 
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VII. APPENDIX: VARIATIONAL WAVE FUNCTION PARAMETERS 

All VAR results presented in the paper have been obtained with light-light correlation functions, F Bq , constructed 
from a linear combination of four gaussians (see eq. (|27[l ). In Table llXl and for the A c baryon with an AL1 inter-quark 
interaction, we present the dependence of the mass and the value of the wave function at the origin on the number of 
gaussians used to build up the correlation function. Finally, we present in Tables Ixl and IXTl the variational parameters 
of the charmed and bottom baryon wave functions for the different quark-quark interactions analyzed in this work. 
The wave functions are easy to handle and can be used to evaluate further observables. 
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AL1 
AL1 X 
AL2 
API 
AP2 
A BD 



AL1 
AL1 X 
AL2 
API 
AP2 
E BD 



AL1 
AL1 X 
AL2 
API 
AP2 
E* BD 



AL1 
AL2 
API 
AP2 
5 BD 



AL1 
AL2 
API 
AP2 
3' BD 



AL1 
AL2 
API 
AP2 
3* BD 



AL1 
AL2 
API 
AP2 
ft BD 



AL1 
AL2 
API 
AP2 
fi* BD 



bi di 



d-2 



Charm 

a 3 b 3 d 3 



Ct4 64 



0.51 
0.57 
1.07 
0.56 
0.72 
1.08 



0.76 
1.51 
0.37 
0.64 
0.77 
0.86 



0.67 0. 
0.68 1. 
0.62 0. 
0.69 0. 
0.79 0. 
0.34 0. 



73 0.73 
11 1.02 
65 0.93 
53 0.62 
72 0.50 
45 0.61 



0.56 0. 
1.31 1. 
0.51 0. 
0.64 0. 
0.78 0. 
0.79 0. 



85 0.76 
49 1.10 
99 0.49 
79 0.68 
46 0.65 
56 0.83 



0.77 1 
0.82 1 
1.13 
0.92 
0.72 
1.28 0. 



12 0.84 
27 0.95 
45 0.94 
92 0.76 
81 0.77 
87 0.74 



0.55 
0.46 
0.71 
0.43 
0.42 
0.58 



-0.12 
0.44 

-0.12 
0.42 
0.43 
0.90 



-0.61 0. 
0.63 0. 
1.00 0. 
0.75 0. 
0.80 0. 
0.99 0. 



62 0.01 
51 0.50 
56 0.50 
51 0.39 
42 0.47 
48 0.67 



0.52 0. 
0.54 0. 
0.51 0. 
0.53 0. 
0.54 0. 
0.76 0. 



63 0.24 
50 0.45 
44 0.26 
38 0.37 
52 0.38 
33 0.71 



0.77 
0.61 
1.17 
0.62 
0.68 
1.08 0. 



50 0.27 

54 0.62 
43 0.68 
49 0.48 

55 0.70 
54 0.65 



0.49 
0.46 
0.67 
0.42 
0.45 
0.75 



0.46 
0.45 
0.44 
0.43 
0.51 
0.87 



0.61 0. 
0.62 0. 
0.92 0. 
0.64 0. 
0.63 0. 
0.93 0. 



50 0.45 
45 0.46 

49 0.70 
48 0.42 

50 0.47 
45 0.68 



0.63 0. 
0.69 0. 
0.86 0. 
0.64 0. 
0.64 0. 
0.75 0. 



53 0.40 

54 0.38 
53 0.05 
46 0.38 
51 0.39 
51 0.58 



0.69 
0.72 
0.90 
0.69 
0.62 
1.22 



0.52 
0.57 
0.48 
0.53 
0.53 
0.65 



0.65 
0.78 
1.10 
0.92 
0.82 



0.86 
0.49 
0.45 
0.71 
0.82 



0.35 0. 
0.82 0. 
0.09 0. 
0.01 0. 
0.45 0. 



73 0.56 
52 0.47 
63 0.47 
69 0.49 
33 0.57 



0.69 0. 
1.06 0. 
0.55 0. 
0.51 0. 
0.77 0. 



50 0.81 
69 0.73 
74 0.51 
87 0.37 
60 0.75 



0.68 1. 
1.65 1 
0.65 
0.54 
2.10 



06 0.62 
12 0.62 
55 0.28 
48 0.43 
86 0.55 



0.53 
0.61 
0.50 
0.61 
0.58 



0.51 
0.46 
0.51 
0.45 
0.78 



0.64 0. 

1.00 0. 

0.66 0. 

0.44 0. 

0.75 0. 



58 0.54 
67 0.47 
57 0.52 
47 0.39 
53 0.61 



0.60 0. 
0.77 0. 
0.60 0. 
0.74 0. 
1.01 0. 



60 0.43 
50 0.39 
60 0.43 
42 0.49 
52 0.62 



0.71 
0.99 
0.69 
0.69 
1.35 



51 0.59 
48 0.41 
50 0.56 
57 0.52 
59 0.61 



0.55 
0.53 
0.49 
0.74 
0.61 



0.47 
0.18 
0.47 
0.38 
0.83 



0.60 0. 
0.71 0. 
0.57 0. 
0.21 0. 
0.69 0. 



55 0.47 
82 0.31 
53 0.44 
57 0.19 
49 0.69 



0.59 0. 
0.71 0. 
0.62 0. 
0.68 0. 
1.09 



54 0.46 
59 0.37 
54 0.38 
57 0.21 
64 0.74 



0.66 
1.36 
0.64 
0.64 
1.48 0. 



53 0.57 
69 0.16 
48 0.50 
47 0.20 
52 0.67 



0.66 
0.71 
0.78 
0.63 
1.02 



0.43 
0.46 
0.14 
0.50 
0.20 



0.10 0. 
0.86 0. 
0.97 0. 
0.58 0. 
0.75 0. 



59 0.58 
75 0.43 

50 0.97 

60 0.60 

51 0.22 



0.92 0. 
0.80 0. 
-0.04 0. 
0.52 0. 
0.68 0. 



56 0.55 
55 0.42 

57 0.80 
67 0.53 

58 0.33 



0.75 0. 
0.83 
0.83 1 
0.76 
1.37 



71 0.65 
69 0.50 
16 1.33 
65 0.64 
90 0.21 



0.74 
0.86 
0.70 
0.67 
0.46 



0.49 
0.17 
0.61 
0.49 
1.10 



0.46 0. 
0.43 0. 
0.40 0. 
0.25 0. 
0.64 0. 



59 0.49 
50 0.47 
61 0.54 
77 0.61 
61 0.77 



0.63 0. 
0.38 0. 
0.75 0. 
0.55 0. 
1.15 0. 



66 0.38 
48 0.53 
60 0.49 
74 0.45 
73 0.75 



0.74 
0.70 
0.72 
0.69 
1.32 



52 0.45 
67 0.62 
54 0.60 
56 0.43 
63 0.77 



V 



0.19 0.19 
0.39 0.39 
0.18 0.18 
0.13 0.13 
0.20 0.20 
0.22 0.22 



0.06 0.06 
0.15 0.15 
0.13 0.13 
0.10 0.10 
0.07 0.07 
0.08 0.08 



0.33 0.33 
0.37 0.37 
0.41 0.41 
0.32 0.32 
0.32 0.32 
0.27 0.27 



0.34 0.15 
0.37 0.25 
0.28 0.25 
0.32 0.20 
0.23 0.18 



0.24 0.22 
0.16 0.04 
0.23 0.20 
0.08 0.08 
0.12 0.11 



0.28 0.28 
0.54 0.40 
0.35 0.36 
0.44 0.50 
0.39 0.36 



0.17 0.17 
0.13 0.13 
0.19 0.19 
0.19 0.19 
0.04 0.04 



0.44 0.44 
0.50 0.50 
0.46 0.46 
0.47 0.47 
0.34 0.34 



TABLE X: Variational parameters (a's are dimensionless, d's have dimensions of fm and b's have dimensions of fm _1 ) of the 
baryon three body wave function (Eqs. I26H and 1271 '). in the charm sector and for different inter-quark interactions. Besides 
for the BD and ALlx interactions, the c parameter (Eq. I28H ) has been set to 200 fm" 1 for the E and E* baryons. In the rest 
of cases, c is set to +00. 
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a 2 



d-2 



Bottom 

a 3 b 3 d 3 



"4 



d 4 



V 



AL1 
AL1 X 
AL2 
API 
AP2 
A BD 



0.60 0.54 
0.48 2.67 
0.45 0.29 
0.56 0.67 
0.96 0.68 
0.78 0.95 



0.58 
1.09 1 
1.21 
0.72 
0.41 
0.99 



71 0.76 
67 1.40 
63 0.91 
53 0.65 
70 0.47 
43 0.76 



0.46 
1.74 
0.96 
0.66 
0.55 
1.22 



75 0.72 
54 1.04 
97 0.62 
80 0.70 

72 0.42 

73 0.79 



0.78 0.84 
0.66 1.30 
1.73 0.88 
0.94 
0.55 
1.24 



0.77 
1.59 
88 0.63 
91 0.78 
43 0.58 
77 0.82 



0.00 
0.09 
-0.08 
0.09 
0.04 
0.04 



0.00 
0.09 
-0.08 
0.09 
0.04 
0.04 



AL1 
AL1 X 
AL2 
API 
AP2 
£ BD 



0.38 
0.49 
0.65 
0.43 
0.40 
0.64 



0.18 
0.44 
0.06 
0.39 
0.46 
0.50 



1.17 0. 
0.68 0. 
0.95 0. 
1.09 0. 
0.71 0. 
1.07 0. 



48 0.44 
47 0.50 
56 0.71 
51 0.44 
45 0.41 
58 0.40 



0.69 
0.59 
0.29 
0.99 
0.71 
1.71 



58 0.33 
54 0.38 
67 0.19 
39 0.41 
50 0.29 
47 0.35 



0.86 0. 
0.65 0. 
1.31 
1.09 0. 
0.75 0. 
1.77 0. 



59 0.74 
46 0.67 
22 0.58 
54 0.55 
37 0.53 
37 0.54 



0.08 
0.21 

-0.08 
0.09 
0.02 

-0.05 



0.08 
0.21 

-0.08 
0.09 
0.02 

-0.05 



AL1 
AL1 X 
AL2 
API 
AP2 
E* BD 



0.44 
0.39 
0.56 
0.39 
0.38 
0.71 



0.45 
0.51 
0.12 
0.51 
0.50 
0.74 



0.63 
0.69 
0.98 
0.70 
0.70 
0.87 



52 0.48 
55 0.40 
45 0.72 
55 0.40 
54 0.39 
31 0.62 



0.58 0. 

0.60 0. 

0.48 0. 

0.59 0. 

0.60 0. 

0.83 0. 



49 0.41 
46 0.41 
55 0.31 
46 0.41 
45 0.41 
40 0.60 



0.68 0. 
0.60 0. 
1.14 0. 
0.60 0. 
0.59 0. 
1.30 0. 



54 0.57 
52 0.45 
32 0.61 
52 0.45 
52 0.45 
44 0.65 



0.23 
0.23 
0.04 
0.22 
0.23 
-0.014 



0.23 
0.23 
0.04 
0.22 
0.23 
-0.014 



AL1 
AL2 
API 
AP2 
3 BD 



0.62 0.61 
0.87 0.40 
0.92 0.63 
0.95 0.74 
0.84 0.85 



0.60 
0.97 
0.47 
0.50 
1.13 



75 0.72 

63 0.58 

64 0.50 
68 0.48 
52 0.86 



0.54 0. 

0.69 0. 

0.71 0. 

0.81 0. 

0.18 0. 



82 0.62 
48 0.56 
74 0.40 
71 0.40 
73 1.06 



0.76 0. 
0.93 
0.68 0. 
0.59 0. 
1.10 1. 



71 0.70 
97 0.73 
46 0.54 
43 0.56 
46 1.12 



0.10 
-0.03 
0.06 
0.00 
0.06 



0.06 
-0.03 

0.04 
-0.05 

0.01 



AL1 
AL2 
API 
AP2 
S' BD 



0.50 0.51 
0.69 0.21 
0.49 0.51 
0.49 0.48 
0.74 0.92 



0.66 
1.12 
0.67 
0.57 
0.90 



57 0.52 

53 0.42 

54 0.52 
52 0.43 
33 0.62 



0.60 
0.58 
0.63 
0.64 
0.89 



60 0.43 
41 0.21 
57 0.46 
53 0.38 
50 0.63 



0.69 0. 
1.29 0. 
0.69 0. 
0.64 0. 
1.26 0. 



50 0.56 

48 0.68 

51 0.59 
47 0.49 

49 0.72 



0.23 
-0.03 

0.23 
-0.04 

0.02 



0.20 
-0.03 

0.20 
-0.03 
-0.05 



AL1 
AL2 
API 
AP2 
H* BD 



0.50 0.51 
0.85 0.06 
0.49 0.52 
0.57 0.48 
0.77 0.86 



0.66 
1.00 
0.68 
0.62 
0.90 



57 0.52 
59 0.50 
51 0.51 
40 0.39 
33 0.63 



0.60 0. 

0.48 0. 

0.65 0. 

0.80 0. 

0.86 0. 



60 0.51 
44 0.26 
53 0.48 
44 0.40 
46 0.63 



0.69 0. 
0.87 0. 
0.70 0. 
0.74 0. 
1.28 0. 



50 0.56 
33 0.67 
52 0.62 
48 0.58 
50 0.71 



0.22 
0.11 
0.23 
0.02 
-0.03 



0.20 
0.11 
0.20 
0.00 
-0.07 



AL1 
AL2 
API 
AP2 
Q BD 



0.60 0.45 

0.62 0.31 

0.72 0.26 

0.63 0.45 

0.82 0.85 



0.58 
0.67 
0.69 
0.78 
0.89 



61 0.58 

65 0.56 

66 0.32 

62 0.58 
43 0.73 



0.49 0. 

0.50 0. 

0.51 0. 

0.49 0. 

0.73 0. 



55 0.50 
59 0.32 
62 0.22 
55 0.50 
54 0.69 



0.72 0. 
0.88 0. 
0.75 0. 
0.72 0. 
1.29 0. 



66 0.66 

65 0.47 
57 0.26 

66 0.66 
56 0.74 



0.13 
0.06 
0.08 
0.13 
-0.07 



0.13 
0.06 
0.08 
0.13 
-0.07 



AL1 
AL2 
API 
AP2 

n* bd 



0.58 0.48 
0.60 0.44 
0.57 0.49 
0.58 0.48 
0.79 0.83 



0.60 0. 
0.62 0. 
0.66 0. 
0.60 0. 
0.84 0. 



59 0.56 
62 0.57 
50 0.50 
59 0.56 
39 0.69 



0.55 0. 

0.53 0. 

0.57 0. 

0.55 0. 

0.75 0. 



57 0.49 
60 0.42 
57 0.48 
57 0.49 
51 0.67 



0.71 0. 
0.77 0. 
0.73 0. 
0.71 0. 
1.32 0. 



61 0.63 

62 0.56 
59 0.62 
61 0.63 
55 0.72 



0.18 
0.13 
0.20 
0.18 
-0.07 



0.18 
0.13 
0.20 
0.18 
-0.07 



TABLE XI: Same as Tabled for the bottom sector. 
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